The Euler-Poincare Equations 
in Geophysical Fluid Dynamics 

Darryl D. Holm 
Theoretical Division and Center for Nonlinear Studies 
Los Alamos National Laboratory, MS B284 
Los Alamos, NM 87545 

dholm@lanl.gov 

Jerrold E. Marsden 
Control and Dynamical Systems 
California Institute of Technology 107-81 
Pasadena, CA 91125 

marsden@cds . caltech.edu 

Tudor S. Ratiu 
Department of Mathematics 
University of California, Santa Cruz, CA 95064 

ratiu@math.ucsc.edu 

To appear in the proceedings volume of the Isaac Newton Institute for Mathematical Sciences 
The Mathematics of Atmosphere and Ocean Dynamics 

April 16, 1998 

Abstract 

Recent theoretical work has developed the Hamilton's-principle analog of Lie-Poisson 
Hamiltonian systems defined on semidirect products. The main theoretical results are 
twofold: 

1. Euler-Poincare equations (the Lagrangian analog of Lie-Poisson Hamiltonian 
equations) are derived for a parameter dependent Lagrangian from a general 
variational principle of Lagrange d'Alembert type in which variations are con- 
strained; 

2. an abstract Kelvin-Noether theorem is derived for such systems. 

By imposing suitable constraints on the variations and by using invariance proper- 
ties of the Lagrangian, as one does for the Euler equations for the rigid body and 
ideal fluids, we cast several standard Eulerian models of geophysical fluid dynamics 
(GFD) at various levels of approximation into Euler-Poincare form and discuss their 
corresponding Kelvin-Noether theorems and potential vorticity conservation laws. The 
various levels of GFD approximation are related by substituting a sequence of velocity 
decompositions and asymptotic expansions into Hamilton's principle for the Euler equa- 
tions of a rotating stratified ideal incompressible fluid. We emphasize that the shared 
properties of this sequence of approximate ideal GFD models follow directly from their 
Euler-Poincare formulations. New modifications of the Euler-Boussinesq equations and 
primitive equations are also proposed in which nonlinear dispersion adaptively filters 
high wavenumbers and thereby enhances stability and regularity without compromising 
either low wavenumber behavior or geophysical balances. 
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1 Introduction 

The Eulerian formulation of the action principle for an ideal fluid casts it into a form that 
is amenable to asymptotic expansions and thereby facilitates the creation of approximate 
theories. This Eulerian action principle is part of the general procedure of the reduction 
theory of Lagrangian systems, including the theory of the Euler-Poincare equations (the 
Lagrangian analog of Lie-Poisson Hamiltonian equations). This setting provides a shared 
structure for many problems in GFD, with several benefits, both immediate (such as a 
systematic approach to hierarchical modeling and versions of Kelvin's theorem for these 
models) and longer term (e.g., structured multisymplectic integration algorithms). 

This paper will be concerned with Euler-Poincare equations arising from a family of 
action principles for a sequence of standard GFD models in purely Eulerian variables at 
various levels of approximation. We use the method of Hamilton's principle asymptotics 
in this setting. In particular, the action principles of these models are related by differ- 
ent levels of truncation of asymptotic expansions and velocity-pressure decompositions in 
Hamilton's principle for the unapproximated Euler equations of rotating stratified ideal 
incompressible fluid dynamics. This sequence of GFD models includes the Euler equa- 
tions themselves, followed by their approximations, namely: Euler-Boussinesq equations 
(EB), primitive equations (PE), Hamiltonian balance equations (HBE), and generalized La- 
grangian mean (GLM) equations. We also relate our approach to the rotating shallow water 
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equations (RSW), semigeostrophic equations (SG), and quasigeostrophic equations (QG). 
Thus, asymptotic expansions and velocity-pressure decompositions of Hamilton's principle 
for the Euler equations describing the motion of a rotating stratified ideal incompressible 
fluid will be used to cast the standard EB, PE, HBE and GLM models of GFD into Euler- 
Poincare form and thereby unify these descriptions and their properties at various levels of 
approximation. See Tables 4.1 and 4.2 for summaries. 

These GFD models have a long history dating back at least to Rossby [1940], Charney 
[1948] and Eliassen [1949], who used them, in their simplest forms (particularly the quasi- 
geostrophic and semigeostrophic approximations), to study structure formation on oceanic 
and atmospheric mesoscales. The history of the efforts to establish the proper equations 
for synoptic motions is summarized by Pedlosky [1987] and Cushman-Roisin [1994]; see also 
Phillips [1963]. One may consult, for example, Salmon [1983, 1985, 1988], Holm, Mars- 
den, Ratiu and Weinstein [1985], Abarbanel, Holm, Marsden, and Ratiu [1986], and Holm 
[1996] for recent applications of the approach of Hamilton's principle asymptotics to derive 
approximate equations in GFD. 

Well before Rossby, Charney, and Eliassen, at the end of the 19th century, Poincare 
[1901] investigated the formulation of the Euler equations for the dynamics of a rigid body 
in Lie algebraic form. Poincare's formulation of the Euler equations for a rigid body carries 
over naturally to the dynamics of ideal continua, as shown by Holm, Marsden and Ratiu 
[1998a], and Poincare's ideas will form the basis of the present study. Some of Poincare's 
other key papers in this area are listed in the bibliography. 

Starting from the action principle for the Euler equations, the present work first expresses 
the various GFD equations in the Euler-Poincare form for continua due to Holm, Marsden 
and Ratiu [1998a] and discusses the properties acquired by casting the GFD equations into 
this form. The main property so obtained is the Kelvin-Noether theorem for the theory. 
This, in turn, leads to conservation of potential vorticity on fluid parcels. Domain-integrated 
energy is also conserved and the relation of the Euler-Poincare equations to the Lie-Poisson 
Hamiltonian formulation of the dynamics is given by a Legendre transformation at the level 
of the Lie algebra of divergenceless vector fields. 

The methods of this paper are based on reduction of variational principles; that is, on 
Lagrangian reduction (see Cendra et al. [1986, 1987] and Marsden and Scheurle [1993a,b]), 
which is also useful for systems with nonholonomic constraints. This has been demonstrated 
in the work of Bloch, Krishnaprasad, Marsden and Murray [1996], who derived the reduced 
Lagrange d'Alembert equations for such nonholonomic systems. Coupled with the methods 
of the present paper, these techniques for handling nonholonomic constraints should also be 
useful for continuum systems. In addition, it seems likely that the techniques of multisym- 
plectic geometry, associated variational integrators, and the multisymplectic reduction will 
be exciting developments for the present setting; see Marsden, Patrick and Shkoller [1997] 
for the beginnings of such a theory. 



Organization of the Paper. In §|| we recall from Holm, Marsden and Ratiu [1998a] 
the abstract Euler-Poincare theorem for Lagrangians depending on parameters along with 
the associated Kelvin-Noether theorem. These theorems play a key role in the rest of our 
analysis. In we discuss their implications for continuum mechanics and then in §|we 
apply them to a sequence of models in geophysical fluid dynamics. We begin in Qand 
§4.2 by recalling the action principles in the Eulerian description for the Euler equations 
and their Eulcr-Boussinesq approximation, respectively. Then we show how these standard 
GFD models satisfy the Euler-Poincare theorem. These sections also introduce the scaling 
regime and small parameters we use in making asymptotic expansions and velocity-pressure 
decompositions that are used in the remaining sections. Next, § [1.3] introduces the hydro- 
static approximation into the Euler-Poincare formulation of the Euler-Boussinesq equations 
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to yield the corresponding formulation of the primitive equations. Later sections cast fur- 
ther approximations of the Euler-Boussinesq equations into the Euler-Poincare formulation, 
starting in §4.4 with the Hamiltonian balance equations and proceeding to the generalized 
Lagrangian mean (GLM) theory for wave, mean flow interaction (WMFI), due to Andrews 
and Mclntyre [1978a,b] in §§. In §| we use the Euler-Poincare theorem, including advected 
parameters, to formulate a new model of ideal GFD called the EBa model that includes 
nonlinear dispersion along with stratification and rotation. The EBa equations modify the 
usual Euler-Boussinesq equations by introducing a length scale, a. The length scale a is 
interpreted physically in the GLM setting as the amplitude of the rapidly fluctuating com- 
ponent of the flow. We derive the Euler-Poincare equations for the EBa model by making 
an asymptotic expansion of the GLM Lagrangian for WMFI in powers of a and the Rossby 
number. Thus, the EBa model is a WMFI turbulence closure model for a rotating stratified 
incompressible fluid. In this model, nonlinear dispersion (parameterized by a) acts to filter 
the high wavenumbers (fc > 1/a) and thereby enhances solution stability and regularity 
without compromising either low wavenumber behavior (fc < 1/a), or geophysical balances. 
We also present the corresponding nonlinear dispersive modification of the primitive equa- 
tions, called the PEa model. The nonlinear dispersive filtering of high wavenumber activity 
in the EBa and PEa models regularizes these equations and thereby makes them good 
candidates for long term numerical integration. 



2 The Euler-Poincare Equations, Semidirect 
Products, and Kelvin's Theorem 

Here we recall from Holm, Marsden and Ratiu [1998a] the statements of the Euler-Poincare 
equations and their associated Kelvin-Noether theorem. In the next section, we will discuss 
these statements in the context of continuum mechanics and then in the following section 
apply them to a sequence of models in geophysical fluid dynamics. Although there are 
several possible permutations of the conventions, we shall state the Euler-Poincare theorem 
for the case of right actions and right invariant Lagrangians, which is appropriate for fluids 
and, in particular, for the GFD situation. 



2.1 The Euler— Poincare Equations and Semidirect Products 

Assumptions and Notation. We shall begin with the abstract framework which will be 
a convenient setting for the several special cases of GFD to follow. 

• Let G be a Lie group and let g be its Lie algebra. We consider a vector space V and 
assume we have a right representation of G on V . The group G then acts in a natural 
way on the right on the dual space V* (the action by g G G on V* is the dual of the 
action by g _1 on V). We denote the action of g on an element v £ V by vg and on 
an element a G V* by ag. In general we use this concatenation notation for group 
actions. Then G also acts by right translation on TG and hence it acts on TG x V*. 
We denote the action of a group element g on a point (vh,a) by (v g , a)g = (vhg, ag). 

• Assume we have a Lagrangian L : TG x V* -^M. that is right G-invariant. 

• For each a Q G V* , define the Lagrangian L ao : TG — > R by L ao {v g ) — L(v g7 a )- Then 
L ao is right invariant under the lift to TG of the right action of G ao on G, where G ao is 
the isotropy group of do (that is, the subgroup of elements of G that leave the element 
oq G V* invariant). 
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• Right G-invariancc of L permits us to define the reduced Lagrangian through the 
equation I : g x V* — ► R by 

l(v g g~ 1 ,ag~ 1 ) = L(v g ,a). 

Conversely, this relation defines for any I : g x V* — > R a right G-invariant function 
L : TG x V* — > R. 

• For a curve p(i) G G, let := .g(t)g(i) _1 and define the curve a(t) — ao.g(t) _1 , 
which is the unique solution of the linear differential equation with time dependent 
coefficients a(t) = —a(t)^(t) with initial condition a(0) = do- 

• Let adj : g — > g be the infinitesimal adjoint operator; that is, the linear map given by 
the Lie algebra bracket: ad^ry) = [£,rj\. Let adj? : g* — > g* be the dual of the linear 
transformation ad ? . 

Theorem 2.1 (Euler Poincare reduction.) The following are equivalent: 
i Hamilton's variational principle 



■L 



L ao (g(t),9(t))dt = (2.1) 



holds, for variations Sg(t) of g(t) vanishing at the endpoints. 
ii g{t) satisfies the Euler-Lagrange equations for L ao on G. 
iii The constrained variational principle 



5 f'l(£(t),a(t))dt = 
Jti 



(2.2) 



holds on g x V* , using variations of the form 

= V ~ adjr? = i] - [€,rj\, 5a = -an, (2.3) 
where rj(t) G vanishes at the endpoints. 
iv The Eulei — Poincare equations hold on g x V* 

dt6£~ tSt + 5a~ 



We refer to Holm, Marsden and Ratiu [1998a] for the proof of this in the abstract setting. 
We shall see some of the features of this result in the concrete setting of continuum mechanics 
shortly. 



Important Notation. Following the notational conventions of Holm, Marsden and Ratiu 
[1998a], we let p v : g — > V be the linear map given by p v {£) = v£ (the right action of £ on 
v e V), and let p* : V* — > g* be its dual. For a £ V* , we write 

p*a = »«s6g', 

which is a bilinear operation in v and a. Continuing to use the concatenation notation for 
Lie algebra actions, the g-action on g* and V* is defined to be minus the dual map of the 
0-action on g and V respectively and is denoted by u£ and a£ for £ G g, a G g* , and a G V*. 
The following is a useful way to write the definition of w o a G g*: for all v G V, a G V* and 
£ G g, we have (note minus sign) 

(voa,0 = (a,vO = - (at,v) . (2.5) 
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The Legendre Transformation. As explained in Marsden and Ratiu [1994], one nor- 
mally thinks of passing from Euler-Poincare equations on a Lie algebra g to Lie-Poisson 
equations on the dual g* by means of the Legendre transformation. In our case, we start 
with a Lagrangian on g x V* and perform a partial Legendre transformation in the variable 
£ only, by writing 

SI 

M =— , h(n,a) = (n,0-l(Z,a). (2.6) 
St; 

Therefore, we have the formulae 

Sh „ I 8£\ /SI 8£\ , , Sh SI /n 

One of the points is that, consistent with the examples, we do not attempt to use the full 
Legendre transformation to make Euler-Poincare equations on g x V correspond to Lie- 
Poisson equations on the dual space 2* x V* . In fact, such attempts will fail because in 
most interesting examples, the full Legendre transform will be degenerate (the heavy top, 
compressible fluids, etc). It is for this reason that we take a partial Legendre transformation. 
In this case, our Euler-Poincare equations on g x V* will correspond to the Lie-Poisson 
equations on g* x V*. We next briefly recall the Hamiltonian setting on g* x V* . 



Lie Poisson Systems on Semidirect Products. Let S = G(§)V be the semidirect 
product Lie group for right actions. Explicitly, the conventions for S are the following: the 
multiplication has the expression 



(ffi)«i)(ff2,V2) = (5152, v 2 +v 1 g 2 ), 



(2.8) 



the identity element is (e,0), and the inverse is given by (g, v) 1 = (g , — vg ). The Lie 
algebra of S is denoted 5 = g (S) V and it has the bracket operation given by 



(£2,^2)] = ([£1,62], «i6 - W2C1). 



(2.9) 



Let s* denote the dual of s. For a right representation of G on V the (+) Lie-Poisson bracket 
of two functions /, k : s* — > R has the expression 



{f,k}+{v,a) = (n 



Sf 5k 
5n 5/j, 



6k Sf Sf 8k \ 
' 8a Sji 8a 8[i / ' 



(2.10) 



where 8f/8u£ g, and Sf/Sa G V are the functional derivatives of /. Using the diamond 
notation (2.5), the corresponding Hamiltonian vector field for h : s* — ■> R is easily seen to 
have the expression 



, . / .„ Sh Sh 

X h (p, a) = - \ ad sh/Sli fx + — o a, a — 



Thus, Hamilton's equations on the dual of a semidirect product are given by 

(1 = {fi, h} = - ad* sh/5fJi fj, - ^oa, 



{a,h\ = - a — 

0/i 



(2.11) 

(2.12) 
(2.13) 



where overdot denotes time derivative. Thus, the partial Legendre transformation (2.6) 
maps the Euler-Poincare equations (2.4), together with the equations a = —a(t)t;(t) for a 
to the Lie-Poisson equations (2.12) and ( 2.13| ). 
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Cautionary Remark. If the vector space V is absent and one has just the equations 



d 61 _ » SI 
dtSj ~ " ad « 6£ 



= -ad*- (2.14) 



for £ e q on a Lie algebra, one speaks of them as the basic Euler-Poincare eq uati ons. 



As explained in Holm, Marsden and Ratiu [1998a], the Euler-Poincare equations (2.4) are 
not the basic Euler-Poincare equations on the larger semidirect product Lie algebra g (S) V * . 
This is a critical difference between the Lie-Poisson and the Euler-Poincare cases. 



Advected Parameters. As we shall see in the examples, and as indicated by the above 
Euler-Poincare reduction theorem, the parameters a £ V* acquire dynamical meaning under 
Lagrangian reduction. For the heavy top, the parameter is the unit vector in the direction 
of gravity, which becomes a dynamical variable in the body representation. For stratified 
incompressible fluids, the parameters are the buoyancy b and volume I? of a fluid element in 
the reference configuration, which in the spatial representation become dynamical variables 
satisfying the passive scalar advection equation and continuity equation, respectively 



2.2 The Kelvin— Noether Theorem 

In this section, we explain a version of the Noether theorem that holds for solutions of 
the Euler-Poincare equations. Our formulation is motivated and designed for continuum 
theories (and hence the name Kelvin-Noether), but it may be also of interest for finite 
dimensional mechanical systems. Of course it is well known that the Kelvin circulation 
theorem for ideal flow is closely related to the Noether theorem applied to continua using 
the particle relabelling symmetry group (see, for example, Arnold [1966]). 



The Kelvin— Noether Quantity. We start with a Lagrangian L ao depending on a pa- 
rameter ao € V* as above. We introduce a manifold C on which G acts (on the right, as 
above) and suppose we have an equivariant map JC : C x V* — > q**. 

In the case of continuum theories, the space C will be a loop space and (/C(c, a), a) for 
c £ C and u G q* will be a circulation. This class of examples also shows why we do not 
want to identify the double dual q** with g. 

Define the Kelvin-Noether quantity J:CxgxV* -tlby 

I(c,Z,a) = (}C(c,a),^(Z,a)). (2.15) 



Theorem 2.2 (Kelvin— Noether) Fixing Co € C, let £(t),a(t) satisfy the Euler-Poincare 
equations and define g(t) to be the solution of g(t) = £{t)g(t) and, say, g(0) — e. Let 
c(t) = cogit)- 1 and I(t) = I(c(t),£(t),a(t)). Then 



j t I(t) = /}C(c(t),a(t)),^oa 



(2.16) 



Proof. First of all, write a(t) = aog(i) 1 and use equivariance to write I(t) as follows: 



IC(c(t),a(t)),^(t),a(t))) =<AV„.»„s. 



9(t) 



The g 1 pulls over to the right side as g (and not g 1 ) because of our conventions of always 
using right representations. We now differentiate the right hand side of this equation. To 
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do so, we use the following well known formula for differentiating the coadjoint action (see 
Marsden and Ratiu [1994], §9.3): 



dt 



W)9{t)] 



9(t), 



where fj, G g*, and £ S Q is given by 

m=9(t)g(t)- 1 . 

Using this and the Euler-Poincare equations, we get 



dt 



— ( /C(c ,a ) 



£(c ,a ), ^ 



^K(*),a(t)) 



git) 



git) 



K.(co,a ), 
/C(c , Oq), 



,*Sl „8l 51 



g(t) 



51 



6a 

IC{co,a )g(ty 1 



o a 



g(t) 



6l_ 

6a. 



lC(c{t),a{t)), 



61 
5 a 



where, in the last steps we used the definitions of the coadjoint action, the Euler-Poincare 
equation and the equivariance of the map /C. H 

Because the advected terms are absent for the basic Euler-Poincare equations, we obtain 
the following. 

Corollary 2.3 For the basic Euler-Poincare equations, the Kelvin quantity I{t), defined 
the same way as above but with I : C x g — > R, is conserved. 



3 The Euler-Poincare Equations in Continuum Mechan- 
ics 

In this section we will apply the Euler-Poincare equations to the case of continuum me- 
chanical systems. 

Vector Fields and Densities. Let T> be a bounded domain in R n with smooth boundary 
DT> (or, more generally, a smooth compact manifold with boundary and given volume form 
or density). We let Diff(2?) denote the diffeomorphism group of T> of an appropriate Sobolev 
class (for example, as in Ebin and Marsden [1970]). If the domain T> is not compact, then 
various decay hypotheses at infinity need to be imposed. Under such conditions, Diff(£>) 
is a smooth infinite dimensional manifold and a topological group relative to the induced 
manifold topology. Right translation is smooth but left translation and inversion are only 
continuous. Thus, Diff(£>) is not literally a Lie group in the naive sense and so the previous 
theory must be applied with care. Nevertheless, if one uses right translations and right 



representations, the Euler-Poincare equations of Theorem 2.1 do make sense, as a direct 
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verification shows. We shall illustrate such computations, by verifying several key facts in 
the proof as we proceed. 

Let X(D) denote the space of vector fields on V of the same differentiability class as 
Diff(£>). Formally this is the right Lie algebra of DifF(D), that is, its standard left Lie 
algebra bracket is minus the usual Lie bracket for vector fields. To distinguish between 
these brackets, we shall reserve in what follows the notation [u, w] for the standard Jacobi- 
Lie bracket of the vector fields u, w G X(D) whereas the notation ad u w := — [u, w] denotes 
the adjoint action of the left Lie algebra on itself. (The sign conventions will also be clear 
in the coordinate expressions.) 

We also let X(T>)* denote the geometric dual space of X(T>), that is, X(D)* := SI 1 (1?) ® 
Den(P), the space of one-form densities on I?. If a <S> m £ O 1 (2?) ® Den(£>), the pairing of 
a ® m with w £ X(T>) is given by 

(a<8>m,w)= / a-wm (3.1) 

where a ■ w denotes the contraction of a one-form with a vector field. For w £ X(T>) and 
a ® m G X(£>)*, the dual of the adjoint representation is defined by (note the minus sign) 

(ad^(a (g) to), u) = — a ■ [w, u] m 
Jv 

and its expression is 

ad^(a (S> to) = (£ w a + (div m w)a) ® to = £ w (a (g) to) , (3.2) 

where div m w is the divergence of w relative to the measure to,, which is related to the Lie 
derivative by ,£ w to = (div m w)m. Hence, if w = w t d/dx l 1 and a — ctidx 1 , the one-form 
factor in the preceding formula for ad^, (a ® to) has the coordinate expression 

/ .deli , dw 1 \ , : / • . Sw-'N , 

the last equality assuming that the divergence is taken relative to the standard measure 
to = rf"x in M™. 

Configurations, Motions and Material Velocities. Throughout the rest of the paper 
we shall use the following conventions and terminology for the standard quantities in contin- 
uum mechanics. Elements of T> representing the material particles of the system are denoted 
by X; their coordinates X A , A = 1, ...,n may thus be regarded as the particle labels. A 
configuration, which we typically denote by 77, is an element of Diff(X>). A motion i] t is 
a path in Diff (V). The Lagrangian or material velocity U(A, t) of the continuum along 
the motion rj t is defined by taking the time derivative of the motion keeping the particle 
labels (the reference particles) X fixed: 

Vt(X), 

x 

the second equality being a convenient shorthand notation of the time derivative holding X 
fixed. 

Consistent with this definition of velocity the tangent space to Diff (2?) at r\ G Diff (2?) 
is given by 

T r) Diff(2?) = {U I( :V^TV\ U V (X) G T„ (X) V}. 



THY t) - d ^ X ^ - d 

upM) Ft 
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Elements of T v DiS(T>) are usually thought of as vector fields on V covering r\. The tangent 
lift of right and left translations on TDiff(£>) by p e Diff(D) have the expressions 

\J r ,i P :=T TI R< p (\J TI )=U ri oip and ip\J v := T v L^(U n ) = Tip o U n . 

Eulerian Velocities. During a motion rj t , the particle labeled by X describes a path in V 
whose points x(X, t) :— r]t(X) are called the Eulerian or spatial points of this path. The 
derivative u(x, t) of this path, keeping the Eulerian point x fixed, is called the Eulerian or 
spatial velocity of the system: 

Thus, the Eulerian velocity u is a time dependent vector field on T>: u t € where 
u t (a;) :— u(x, t). We also have the fundamental relationship 

U t = Ut O T) t , 

where U t (X) := XJ(X,t). 

The representation space V* of Diff(P) in continuum mechanics is often some subspace 
of 1(T>) (§} Den(X>), the tensor field densities on V and the representation is given by pull 
back. It is thus a right representation of Diff(2?) on 1(V) (g)Den(X>). The right action of the 
Lie algebra 3£(D) on V* is given by au := £ u a, the Lie derivative of the tensor field density 
a along the vector field u. 

The Lagrangian. The Lagrangian of a continuum mechanical system is a function L : 
TDiff(P) xV'-tl which is right invariant relative to the tangent lift of right translation 
of Diff(P) on itself and pull back on the tensor field densities. 

Thus, the Lagrangian L induces a reduced Lagrangian I : X(T>) xF-tl defined by 

^(u, a) = L(u o 77, rj*a), 

where u e X(T>) and a e V* C X(2?) <8> Den(P), and where rj*a denotes the pull back of a by 
the diffeomorphism r\ and u is the Eulerian velocity. The evolution of a is given by solving 
the equation 

a = —£ u a. 

The solution of this equation, given the initial condition do, is a(t) — tpt*ao, where the lower 
star denotes the push forward operation and ip t is the flow of u. 

Advected Eulerian Quantities. These are defined in continuum mechanics to be those 
variables which are Lie transported by the flow of the Eulerian velocity field. Using this 
standard terminology, the above equation states that the tensor field density a (which may 
include mass density and other Eulerian quantities) is advected. 

As we have discussed, in many examples, V* C 1(1?) cg>Den(2?). On a general manifold, 
tensors of a given type have natural duals. For example, symmetric covariant tensors are 
dual to symmetric contravariant tensor densities, the pairing being given by the integration 
of the natural contraction of these tensors. Likewise, fc-forms are naturally dual to (n — k) 
forms, the pairing being given by taking the integral of their wedge product. 



u(x,t) :=U(X,t) := - 
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The operation o between elements of V and V* producing an element of £(X>)* intro- 
duced in equation (2.5) becomes 



(voa,w)— / a-£ v/ v = — / v-f^a 
Jv Jt> 



(3.3) 



where v ■ £ w a denotes the contraction, as described above, of elements of V and elements 
of V* . (These operations do not depend on a Riemannian structure.) 

For a path rjt <E D iff (P ) let u(x, t) be its Eulerian velocity and consider, as in the 
hypotheses of Theorem |2.l| the curve a(t) with initial condition ao given by the equation 



a + £ u a = 0. 



(3.4) 



Let L ao (U) := L(U, ao). We can now state Theorem 2A in this particular, but very useful, 
setting. 

Theorem 3.1 (Euler Poincare reduction for continua.) For apathy mDiff(P) with 
Lagrangian velocity U and Eulerian velocity u, the following are equivalent: 

i Hamilton's variational principle 



L(X,U t (X),a Q (X)) dt = Q 

holds, for variations Srj t vanishing at the endpoints. 
ii r) t satisfies the Euler-Lagrange equations for L ao on Diff (T>) . 
iii The constrained variational principle in Eulerian coordinates 

l(u,a) dt = 



(3.5) 



holds on X(T>) x V* , using variations of the form 

. dw dw . . 

5u = — - ad u w = — + [u,wj, 

where w t = dr]t o rj^ vanishes at the endpoints. 
iv The Eulei — Poincare equations for continua 



5a 



(3.6) 



(3.7) 



61 61 



d 81 „ 51 51 

■7777- = -ad u — + — oa = -£uy~ + T oa > 
at ou 5u da ou da 



(3.8) 



hold, where the o operation given by (S.i) needs to be determined on a case by case 
basis, depending on the nature of the tensor a. (Remember that 51 /5u is a one-form 
density.) 



Remarks. 

1. The following string of equalities shows directly that iii is equivalent to iv: 







51 51 
l(u, a)dt = I ( — • 5u + — ■ 5a ) dt 
' 1 ou 5a 



pt2 


'81 


( dw 


It! 


5u 


' \~~dt ~ 


/"*2 

/ w • 

'tl 


d 61 


~~ dt6u ' 



ad„ w 



51 
8a 
51 51 
8u 5a 



dt 



dt 



(3.9) 
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Similarly, one can deduce the form (3.7) of the variations in the constrained variational 
principle ( |3.6[ ) by a direct calculation, as follows. One writes the basic relation between 
the spatial and material velocities, namely u(x,t) = r](ri^ 1 (x),t). One then takes 
the variation of this equation with respect to rj and uses the definition w(x, t) = 
5rj((7] t T 1 (x),t) together with a calculation of its time derivative. Of course, one can 
also do this calculation using the inverse map rj^ 1 instead of n as the basic variable, 
see Holm, Marsden, and Ratiu [1986], Holm [1996]. 

The preceding sort of calculation for Su in fluid mechanics and the interpretation of 
this restriction on the form of the variations as the so-called Lin constraints is due to 



Brctherton [1970]. The variational form (3.7) for the 'basic' Euler-Poincare equations 



(i.e., without the advected parameters a) is due to Marsden and Scheurle [1993a] and 
Bloch, Krishnaprasad, Marsden and Ratiu [1996]. For the full Euler-Poincare case, 
this form is due to Holm, Marsden and Ratiu [1998a] and for the general Lagrangian 
reduction case, to Cendra, Holm, Marsden and Ratiu [1997] and Cendra, Marsden and 
Ratiu [1997]. These ideas were used for Maxwell- Vlasov plasmas by Cendra, Holm, 
Hoyle and Marsden [1997]. 

4. The coordinate expressions for (51 / Sa)oa required to complete the equations of motion 
for GFD models are given in the next section for ao(X) in three dimensions. 

5. As with the general theory, variations of the action in the advected tensor quantities 
contribute to the equations of motion which follow from Hamilton's principle. At the 
level of the action I for the Euler-Poincare equations, the Legendre transform in the 
variable u alone is often nonsingular, and when it is, it produces the Hamiltonian 
formulation of Eulerian fluid motions with a Lie-Poisson bracket defined on the dual 
of the semidirect product algebra of vector fields acting amongst themselves by Lie 
bracket and on tensor fields and differential forms by the Lie derivative. This is a 
special instance of the more general facts for Lie algebras that was discussed earlier. 

The Inverse Map and the Tensor Fields a for Fluids. In the case of fluids in the 
Lagrangian picture, the flow of the fluid is a diffeomorphism which takes a fluid parcel along 
a path from its initial position X, in a "reference configuration" to its current position 
x in the "container", i.e., in the Eulerian domain of flow. As we have described, this 
forward map is denoted by rj : X i— > x. The inverse map ry _1 : x i— > X provides the 
map assigning the Lagrangian labels to a given spatial point. Interpreted as passive scalars, 
these Lagrangian labels are simply advected with the fluid flow, X = 0. In the Lagrangian 
picture, a tensor density in the reference configuration ao(X) (satisfying clq(X) = 0) consists 
of invariant tensor functions of the initial reference positions and their differentials. These 
tensor functions are parameters of the initial fluid reference configuration (e.g., the initial 
density distribution, which is an invariant n-form). 
When viewed in the Eulerian picture as 

a t (x) := (r, t *a )(x) = (% 1 *a )(x), 

i.e., 

ao(X) := (r)* t a t )(X) = {^a ){X), 

the time invariant tensor density ao(X) in the reference configuration acquires advective 
dynamics in the Eulerian picture, namely 

do(X) = (§- t + £») a(x,i)=0, 
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where £ u denotes Lie derivative with respect to the Eulerian velocity field u(x,t). This 
relation results directly from the well known Lie derivative formula for tensor fields. (See, 
for example, Abraham, Marsden and Ratiu [1988].) 

Mapping the time invariant quantity ao (a tensor density function of X) to the time 
varying quantity a t (a tensor density function of x) as explained above is a special case of 
the way we advect quantities in V* in the general theory. Specifically, we can view this 
advection of at as being the fluid analogue of the advection of the unit vector along the 
direction of gravity (a spatially fixed quantity) by means of the body rotation vector in the 
heavy top example. 

Consistent with the fact that the heavy top is a left invariant system while continuum 
theories are right invariant, the advected tensor density a t is a spatial quantity, while the 
advected direction of gravity is a body quantity. If we were to take the inverse map rj^ 1 as 
the basic group variable, rather than the map 77, then continuum theories would also become 
left invariant. 



The Continuity Equation for the Mass Density. We will need to make an additional 
assumption on our continuum theory. Namely, we assume that amongst the tensor densities 
being advected, there is a special one, namely the mass density. This of course is a tensor 
density that occurs in all continuum theories. We denote this density by pd n x and it is 
advected according to the standard principles discussed above. Thus, p satisfies the usual 
continuity equation: 

d 

— p+div(pu) = 0. 

In the Lagrangian picture we have pd n x = po(X)d n X, where pa(X) is the mass density in 
the reference configuration. It will also be convenient in the continuum examples below to 
define Lagrangian mass coordinates t{X) satisfying pd n x = d n £ with I = 0. (When using 
Lagrangian mass coordinates, we shall denote the density p as D.) 



The Kelvin Noether Theorem. Let 

C := {7 : S 1 — > T> I 7 continuous} 

be the space of continuous loops in D and let the group Diff(2?) act on C on the left by 
(77, 7) £ Diff(2?) x C 1— > 777 £ C, where 777 = 77 o 7. 

Next we shall define the circulation map JC : C x V* — > X(T>)** . Given a one form 
density a 6 X* we can form a one form (no longer a density) by dividing it by the mass 
density p; we denote the result just by a/ p. We let JC then be defined by 

(£(7,a),a)= (3.10) 

The expression in this definition is called the circulation of the one-form a/p around the 
loop 7. 

This map is equivariant in the sense that 

(JC(rj o 7, ri*a),ri*a} = (£(7, a), a) 

for any 77 € Diff (T>). This is proved using the definitions and the change of variables formula. 

Given the Lagrangian I : X(V) x V* — > R, the Kelvin-Noether quantity is given by ( 2.15| ) 
which in this case becomes 

/(7t,u,a)= / ~. (3.11) 
H P Su 
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With these definitions of K, and /, the statement of Theorem 2.2 becomes the classical Kelvin 
circulation theorem. 



Theorem 3.2 (Kelvin Circulation Theorem.) Assume that u(x,t) satisfies the Euler- 
Poincare equations for continua (3.8): 



d ( 51 



dt \Su 



SI 



51 



and a satisfies 



da 

m 



Su 



+ £ u a = 0. 



So, 



Let rjt be the flow of the Eulerian velocity field u, that is. Ut — (drjt/dt) o rj t 1 . 
It ■= Vt ° 7o and I(t) := I(j t ,u t , a t ). Then 



d . , I 151 
—I(t) = f -j-oa 

dt J lt p 5a 



Define 



(3.12) 



In this statement, we use a subscript t to emphasise that the operations are done at a 
particular t and to avoid having to write the other arguments, as in a t (x) = a(x, t); we omit 
the arguments from the notation when convenient. Due to the importance of this theorem 
we shall give here a separate proof tailored for the case of continuum mechanical systems. 



Proof. First we change variables in the expression for I(t): 



I(t) 



Vt 



1 51 
, ltPt Su JJ0 

Next, we use the Lie derivative formula, namely 



1 51 
Pt Su 



Po 



-Vt 



51 
~6u~ 



Vt 



9 _i_ e 
—a t + £ u a t 

at 



applied to a one-form density at- This formula gives 

d r i 



dt m 



dt ./ 70 po Vt 
l_d_ f , 

70 Po dt \ l 

1 , 
— Vt 

70 P0 



5l_ 
5u 



d_ 

dt 



51 
5u 
5l_ 
5u 



OU 



By the Euler-Poincare equations (|3.8|), this becomes 



dt 



I(t) 



1 

Po 



-Vt 



61 

5a. 



o a 



1 

it Pt 



5a 



o a 



again by the change of variables formula. I 

Corollary 3.3 Since the last expression holds for every loop jt, we may write it as 



d 



dt 



1 51 151 



p 5u p 5a 



(3.13) 
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4 Applications of the Euler— Poincare Theorem in GFD 



Variational Formulae in Three Dimensions. We compute explicit formulae for the 
variations 5a in the cases that the set of tensors a is drawn from a set of scalar fields and 
densities on R 3 . We shall denote this symbolically by writing 



e {b,Dd 3 x}. 



(4.1) 



We have seen that invariance of the set a in the Lagrangian picture under the dynamics of 
u implies in the Eulerian picture that 



5 + * 



0. 



where £ u denotes Lie derivative with respect to the velocity vector field u. Hence, for a 
fluid dynamical action & = J dt l(u; b, D), the advected variables {b, D} satisfy the following 
Lie-derivative relations, 



d \ db 
- + A,j6 = 0, or ^=-u.V6, 

— +£ u )Dd 3 x = 0, or _ = -V-(Du) 



(4.2) 
(4.3) 



In fluid dynamical applications, the advected Eulerian variables b and D represent the 
buoyancy b (or specific entropy, for the compres sibl e case) and vo lume element (or mass 
density) D, respectively. According to Theorem 3J, equation (3.6), the variations of the 
tensor functions a at fixed x and t are also given by Lie derivatives, namely 5a = — £ VJ a, or 



=£ w b = — w • V b , 
£ v , (D d 3 x) = - V • (Dw) d 3 x . 



5D d 3 x = 

Hence, Hamilton's principle with this dependence yields 



(4.4) 







f J dt l{u;b,D) 



51 
SI 



dt 
dt 
dt w • 
dt w 



5u 



61 
Sb 



6b - 



51 
6D 



5D 



/<9 W , \ 51 _, SI i ^_ 
{^- ad ^)-5b W - Vb -5D^- (Dw) 



dt Su u Su 



51 51 
— Vb + D V — 



d \6l 51 , „ _ SI' 
dt +£ »)6^ + Tb Vb - DV 5D 



(4.5) 



where we have consistently dropped boundary terms arising from integrations by parts, by 
invoking natural boundary conditions. Specifically, we impose n • w = on the boundary, 
where n is the boundary's outward unit normal vector. 



The Euler-Poincare equatio ns f or continua (3.8) may now be summarized for advected 
Eulerian variables a in th e set (LI). We adopt the notational convention of the circulation 
map K, in equation ( 3.10 ) that a one form density can be made into a one form (no longer 
a density) by dividing it by the mass density D and use the Lie-derivative relation for the 
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continuity equation (d/dt + £ u )Dd 3 x = 0. Then, the Euclidean components of the Euler- 



Poincare equations for continua in equation (4.5) are expressed in Kelvin theorem form 
(3.13) with a slight abuse of notation as 



1 51 



1 Sl„, 



dx 



y(-) 



dx = 0, 



(4.6) 



in which the variational derivatives of the Lagrangian I are to be computed according to 
the usual physical conventions, i.e., as Frechet derivatives. Formula (4.6) is the Kelvin- 
Noether form of the equation of motion for ideal continua. Hence, wc have the explicit 
Kelvin theorem expression, cf. equations (3.11) and (3.12), 



d_ 
dt 



7t(u) 



1 51 , 
D 5u 



7t(u) 



1 si„, , 
DTb Vb - dx > 



(4.7) 



where the curve 7t(u) moves with the fluid velocity u. Then, by Stokes' theorem, the Euler 
equations generate circulation of (77 Jj) whenever V6 and V(-jylr) are n °t collinear. The 
corresponding conservation of potential vorticity q on fluid parcels is given by 



Oq 
di 



— — h u • Vq = , where q 



1 ™ , / 1 <« 
— V6-curl — — 

D \D5u 



(4. 



Remark on Lagrangian Mass Coordinates. An alternative way to treat Hamilton's 
principle for an action & = J dt l(u;b,D) is to perform variations at fixed x and t of 
the inverse maps x — > i, described by the Lagrangian mass coordinate functions £ A (x,t), 
A = 1,2, ... ,n, which determine u, b and D by the formulae (in which one sums on repeated 
indices) 



d£_ A 
~dt 



ftp A 

— = -u i Df, b = b(£ A ), Df = ^, D = dct(Df) 



(4.9) 



As discussed above, the relation of mass coordinates £ to the usual Lagrangian coordinates 
X is given by a simple change of variables in the fluid reference configuration to make 
po(X)d n X = d n l. Variation of an action of the form 6 = J dt Z(u, b, D) with respect to i A 
with p imposing volume preservation then yields (Holm, Marsden, and Ratiu [1986], Holm 
[1996]), 



SS = 



dt / d 



[d{D- 1 )a 

1 51 
+ DTb> 



5t 



d 1 51 
dtD&u 1 

\6DJa 



1 51 



Sp(D-l)} 



(4.10) 



where d/dt — d/dt + u • V is the material derivative and we again invoke natural boundary 
conditions (n • u = on the boundary) when i ntegra ting by parts. Hence, the vanishing of 
the coefficient of 5l A in the variational formula ( 4.10 ) r ecove rs t he E uler-Poincare equations 
for continua (|3.8| ) expressed in Kelvin theorem form ( 3.13 ) or (4/7), by stationarity of the 
action 6 = J dt l(u;b,D) with respect to variations of the Lagrangian mass coordinates 
l A (x, t). In vector notation, these equations are 



d 1 51 I 51 
dtDSu + DM 

or, in three dimensions. 



Vu 3 



1 51. 



-V6 



51 

V ^^°' 



d_(}_5l_ 
dt\D5u 



— u x curl 



1 51 
D5u 



51 \ 

Id) 



1 51 
DSb 



Vb = 0. 



(4.11) 



(4.12) 
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In writing the last equation, we have used the fundamental vector identity of fluid dynamics, 



(b- V)a- 



b x (V x a) + V(b ■ a) 



(4.13) 



for any three dimensional vectors a and b with, in this case, a = (t>M:) and b = u. Taking 
the curl of equation 



.12 ) and using advection of the buoyancy b yields conservation of 
potential vorticity on fluid parcels as given in equation (1.8). For incompressible flows D = 1 
in equation (|4.S|). The Euclidean component formulae ( 4.11 ), ( 4.12| ) and (IS) are especially 
convenient for direct calculations of motion equations in geophysical fluid dynamics, to 
which we turn our attention next. 



4.1 Euler's Equations for a Rotating Stratified Ideal Incompress- 
ible Fluid 

The Lagrangian. In the Eulerian velocity representation, we consider Hamilton's princi- 
ple for fluid motion in an three dimensional domain with action functional S = J dtl and 
Lagrangian Z(u, b, D) given by 



I = 



J d 3 x p D(l + b)fyu\ 2 + u-R(x)- 3 z) -p(D-l), 



(4.14) 



where ptot = PoD(l + b) is the total mass density, po is a dimensional constant and R is a 
given function of x. This Lagrangian produces the following variations at fixed x and t 



D~5u 
5l_ 
SD 



po(l + 6)(u + R), 



5I_ 

5l> 



PoD 



Po 



(l + fe)(i|u| 2 + u.R 



gz 



l 

2 
P, 



u|" + u ■ R — gz 



51 
5p 



(D-l). 



(4.15) 



Hence, from the Euclidean component formula (4.11) for Hamilton principles of this type 
and the fundamental vector identity ( 4.13| ), we find the motion equation for an Euler fluid 
in three dimensions, 



du 



u x curl R + gz ■ 



1 



Po(l + b) 



Vp = 0, 



(4.16) 



where curlR = 2J~2(x) is the Coriolis parameter (i.e., twice the local angular rotation fre- 
quency) . In writing this equation, we have used advection of buoyancy, 



db 
di 



u • V6 = 0, 



from equation (4.2). 



The Kelvin— Noether Theorem. From equation (4.7), the Kelvin-Noether circulation 
theorem corresponding to the motion equation (4.16) for an ideal incompressible stratified 
fluid in three dimensions is, 



dt 



— f (u + R) • dx 



7t(u) 



1 



t (u) Po(l + b) 



Vp ■ dx 



(4.17) 



where the curve 7t(u) moves with the fluid velocity u. By Stokes' theorem, the Euler 
equations generate circulation of (u + R) around 74 (u) whenever the gradients of bouyancy 
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and pressure are not collinear. Using advection of buoyancy 6, one finds conservation of 
potential vorticity q^ui on fluid parcels, cf. equation (4.8), 



dq 



Eul 



dt 



Vq 



Eul 



, where qEui = V6 ■ curl (u + R) . 



(4.18) 



The constraint D — 1 (volume preservation) is imposed by varying p in Hamilton's prin- 
ciple, according to equation ( 4.15 ). Incompressibility then follows from substituting D = 1 
into the Lie-derivative re lation ([1.3]) for D, which gives Vu = 0. Upon taking the divergence 
of the motion equation (4.16) and requiring incompressibility to be preserved in time, one 
finds an elliptic equation for the pressure p with a Neumann boundary condition obtained 
from the normal component of the motion equation (4.16) evaluated on the boundary. 



4.2 Euler-Boussinesq Equations 

The Lagrangian. The Lagrangian ( [4.14 ) for the Euler fluid motion nondimensionalizes as 
follows, in terms of units of L for horizontal distance, Bq for vertical depth, Uq for horizontal 
velocity, BqUq/L for vertical velocity, fo for Coriolis parameter, po for density and pofoLUo 
for pressure: 



d 3 x D(l + b) ( -u 3 • v 3 + u • R(x) 



p(D-l) 



(4.19) 



Here we take R(x) to be horizontal and independent of the vertical coordinate, so curl R = 
f(x.)z. In this nondimensional notation, three-dimensional vectors and gradient operators 
have subscript 3, while horizontal vectors and gradient operators are left unadorned. Thus, 
we denote, in three dimensional Euclidean space, 



x 3 
u 3 = 



= (x,y,z), x 
(u, v, w), u = 



^_ d d d 

dx' dy 1 dz 



(u,v,0), 

V - (— — o 

\dx' dy' 



d 
dt 



d_ 

dt 



u 3 • V 3 



dt 



u-V 



w 



dz 



(4.20) 



and z is the unit vector in the vertical z direction. For notational convenience, also we 
define two nondimensional velocities 



u 3 = (u, v,w), and v 3 = (u, v, <j 2 w) , 



as well as three nondimensional parameters, 



e = 



U 
foL' 



gB 



(4.21) 



(4.22) 



corresponding respectively to Rossby number, aspect ratio and (squared) rotational 
Froude number. Typically, the Rossby number and the aspect ratio are small, e,er -c 1, 
while the rotational Froude number is of order unity in atmospheric and oceanic dynamics 
for L at synoptic scales and larger. The no ndim ensional Euler fluid equations cor- 
responding to the Lagrangian I in equation (4.19) are obtained from the Euler-Poincare 
equations (4.11) or (4.12) with u — * u 3 and V — > V 3 . Namely, 



dv 3 i r> 

e u x curlR 

dt 



(1 + b) 



V 3 p = 0. 



(4.23) 



18 



Clearly, the leading order balances in these equations are hydrostatic in the vertical and 
geostrophic in the horizontal direction. In this notation, the Kelvin-Noether circulation 
theorem ((0) for the Euler fluid becomes 



— I (ev 3 + R) • dx 3 = — <f ,} ,, V 3 p • dx 3 , 

™/t,(u 3 ) / 7t (u 3 ) Po(1 + &) 



(4.24) 



Likewise, conservation of nondimensional potential vorticity gEui on fluid parcels is given 
by, cf. equation ( 4.18 ), 



dt 



u 3 • V 3 g Eu i = , where g Eu i = V 3 6 ■ V 3 x (ev 3 + R) 



(4.25) 



Hamilton's Principle Asymptotics. For sufficiently small buoyancy, b = o(e), we define 

z , . , b 



P =P + 



and b' 



and expand the Lagrangian (4.19) in powers of e as 



EB 



J dt J d 3 x d(Ui 3 ■ v 3 + u- R(x) -b'zj -p'(D - l)+o(e). 



(4.26) 



Upon dropping the order o(e) term in the Lagrangian Zeb the corresponding Eulcr-Poincare 
equation gives the Euler- Boussinesq equation for fluid motion in three dimensions, 
namely, 



e ^Xi _ u x curl R + b 'z + V 3 p ' 
dt 



0. 



or, in horizontal and vertical components, with curlR = f(x)z 
du 



e — + fz x u + Vp' = 0, 
at 



2 dw , , dp' 

+ 6' + ^-=0, 



dt 



dz 



(4.27) 



(4.28) 



where 



^t- = and 

dt 



dw 



V 3 -u 3 = V- u+ — = 0. 

dz 



Even for order O(e) buoyancy, the leading order bal ance s are still hydrostatic in the ver- 
tical, and geostrophic in the horizontal. Equations ( 4.28 ) describe the motion of an ideal 
incompressible stratified fluid relative to a stable hydrostatic equlibrium in which the den- 
sity is taken to be constant except in the buoyant force. See, for example, Phillips [1969] 
for a derivation of this approximate system via direct asymptotic expansions of the Euler 
equations. 



The Kelvin— Noether Theorem. Th e Kelvin-Noether circulation theorem (4.7) for the 
Euler-Boussinesq motion equation (4.27) for an ideal incompressible stratified fluid in three 
dimensions is, 



d_ 

dt 



(ev 3 + R) • dx = — (p b 'dz , 

7t("3) ^7t(u3) 



(4.29) 



where the curv e 7t(u g_) m oves with the fluid velocity u 3 . (The two Kelvin theorems in 
equations ( 4.24 ) and ( 4.29 ) differ in their right hand sides.) By Stokes' theorem, the Eulcr- 
Boussinesq equations generate circulation of ev 3 + R around 74(113) whenever the gradient 
of bouyancy is not vertical. Conservation of potential vorticity g E B on fluid parcels for the 
Euler-Boussinesq equations is given by 

+ u 3 • V 3 5eb — , where q EB = V 3 & ' • V 3 x (ev 3 + R) . (4.30) 
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4.3 Primitive Equations 

The Lagrangian. The primitive equations (PE) arise from the Euler Boussinesq equa- 
tions, upon imposing the approximation of hydrostatic pressure balance. Setting the aspect 
ratio parameter a to zero in the Lagrangian ^eb in equation ( [f.26| ) (see equations (4.21) and 
(4.22)), provides the Lagrangian for the nondimensional primitive equations (PE), 



^pe = / dt 



(fx 



-D|u| 2 + L>u- R(x) - Db'z-p'(D - 1) 



(4.31) 



The Euler-Poincare equations for now produce the PE; namely, equations ( 4.2q ) with 
a = 0, 



du , , dp' 

e— +/zxu + V P ' = 0,ft'+f - =0, 
at oz 



(4.32) 



where 



aV_ 
dt 



and V3 • u 3 = V • u 



dw 



= 0. 



Thus, from the viewpoint of Hamilton's principle, imposition of hydrostatic balance corre- 
sponds to ignoring the kinetic energy of vertical motion by setting a = in the nondimen- 
sional EB Lagrangian ( 1.26| ). 



The Kelvin— Noether Theorem. The Kelvin-Noether circulation theorem for the prim- 



itive equations is obtained from equation (4.29) for the Euler-Boussinesq equations simply 
by setting a = 0. Namely, 



^r- (f (eu + R) • dx = - (f b'dz 



(4.33) 



where the curve 74(113) moves with the fluid velocity 113. By Stokes' theorem, the primitive 
equations generate circulation of eu + R around 74(113) whenever the gradient of bouyancy 
is not vertical. The conservation of potential vorticity on fluid parcels for the primitive 
equations is given by, cf. equation (4.5 ), 



dq 



PE 



dt 



U3 ■ V 3 <7pe = , where <?pe = V36 ' • V 3 x (eu + R) 



(4.34) 



Remark. In the limit, e — > 0, Hamilton's principle for either Ieb, or Zpe gives, 

fz x u + b'z + V 3 p' = 0, (4.35) 

which encodes the leading order hydrostatic and geostrophic equilibrium balances. These 
balances form the basis for further approximations for near-geostrophic, hydrostatic flow. 



4.4 Hamiltonian Balance Equations 

Balanced Fluid Motions. A fluid motion equation is said to be balanced, if specification 
of the fluid's stratified buoyancy and divergenceless velocity determines its pressure through 
the solution of an equation which does not contain partial time-derivatives among its highest 
derivatives. This definition of balance makes pressure a diagnostic variable (as opposed to 
the dynamic, or prognostic variables such as the horizontal velocity components). The Euler 
equations (4.23) and the Euler-Boussinesq equations (4.27) for the incompressible motion 
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of a rotating continuously stratified fluid are balanced in this sense, because the pressure 
in these cases is determined diagnostically from the buoyancy and velocity of the fluid by 
solving a Neumann problem. However, the hydrostatic approximation of this motion by the 
primitive equations (PE) is not balanced, because the Poisson equation for the pressure in PE 
involves the time-derivative of the horizontal velocity divergence, which alters the character 
of the Euler system from which PE is derived and may lead to rapid time dependence, 
as discussed in Browning et al. [1990]. Balanced approximations which eliminate this 
potentially rapid time dependence have been sought and found, usually by using asymptotic 
expansions of the solutions of the PE in powers of the small Rossby number, e«l, after 
decomposing the horizontal velocity u into order O(l) rotational and order 0(e) divergent 
components, as u = z x Vip + eVx, where tp and \ are the stream function and velocity 
potential, respectively, for the horizontal motion. (This is just the Helmholtz decomposition 
with relative weight e.) 

Balance equations (BE) are reviewed in the classic paper of McWilliams and Gent [1980]. 
Succeeding investigations have concerned the well-posedness and other features of various 
BE models describing continuously stratified oceanic and atmospheric motions. For example, 
consistent initial boundary value problems and regimes of validity for BE are determined in 
Gent and McWilliams [1983a, b]. In other papers by these authors and their collaborators 
listed in the bibliography, balanced models in isentropic coordinates are derived, methods 
for the numerical solution of BE are developed, and the applications of BE to problems of 
vortex motion on a /?-plane and wind-driven ocean circulation are discussed. In studies of 
continuously stratified incompressible fluids, solutions of balance equations that retain terms 
of order O(l) and order O(e) in a Rossby number expansion of the PE solutions have been 
found to compare remarkably well with numerical simulations of the PE; see Allen, Barth, 
and Newberger [1990a, b] and Barth et al. [1990]. Discussions of the relation between BE 
and semigeostrophy have also recently appeared, see, e.g., Gent, McWilliams and Snyder 
[1994] and Xu [1994]. 

Conservation of Energy and Potential Vorticity. One recurring issue in the liter- 
ature is that, when truncated at order 0(e) in the Rossby number expansion, the BE for 
continuously stratified fluids conserve energy (Lorenz [I960]), but do not conserve potential 
vorticity on fluid parcels. Recently, Allen [1991] found a set of BE for continuously stratified 
fluids that retains additional terms of order 0(e 2 ) and does conserve potential vorticity on 
fluid parcels. Allen calls these balance equations "BEM equations", because they are based 
on momentum equations, rather than on the equation for vertical vorticity, as for the stan- 
dard BE. An advantage of the momentum formulation of BEM over the vorticity formulation 
of the original BE is that boundary conditions are more naturally imposed on the fluid's 
velocity than on its vorticity. Holm [1996] derives Hamiltonian balance equations (HBE) 
in the momentum formulation by using the e-weightcd Helmholtz decomposition for u and 
expanding Hamilton's principle (HP) for the PE in powers of the Rossby number, f«l. 
This expansion is truncated at order O(e), then all terms are retained that result from tak- 
ing variations. As we have seen, an asymptotic expansion of HP for the Eulcr-Boussinesq 
(EB) equations which govern rotating stratified incompressible inviscid fluid flow has two 
small dimensionless parameters: the aspect ratio of the shallow domain, a, and the Rossby 
number, e. Setting a equal to zero in this expansion yields HP for PE. Setting e also equal 
to zero yields HP for equilibrium solutions in both geostrophic and hydrostatic balance. 
Setting a = 0, substituting the e- weighted Helmholtz decomposition for u and truncating 
the resulting asymptotic expansion in e of the HP for the EB equations, yields HP for a set of 
nearly-geostrophic Hamiltonian balance equations (HBE). The resulting HBE are equivalent 
to the BEM equations in Allen [1991]. 
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The Lagrangian. The Lagrangian for the HBE model is given in Holm [1996], cf. equa- 
tion (|4.31) for the PE action, 



6hbe = J dt J d 3 x Du-TL(x)- Dbz-p(D -1) + e^\u-eu D \ 2 , (4.36) 

where the horizontal fluid velocity is taken in balance equation form as u = u R + eup = 
z x + eVx- The corresponding Euler-Poincare equations give the dynamics of the HBE 
model 

e^u R + e 2 u Rj Vu 3 D + fz x u + Vp = 0, 

dp 2 du D 
b+ — +e u R - = 0, 

oz oz 

, db d , , db 

with — = — b + u • Vo + e w— = 0, 

dt at oz 

3w 

and V-u + e— = 0. (4.37) 

oz 

Here the notation is the same as for the PE, except that w — > e w for HBE. 

Dropping all terms of order 0(e 2 ) from the HBE model equations fl4.37| ) recovers the BE 
discussed in Gent and McWilliams [1983a, b]. Retaining these order 0(e 2 ) terms restores 
the conservation laws due to symmetries of HP at the truncation order O(e). As explained 
in Holm [1996], the resulting HBE model has the same order 0(e) accuracy as the BE, since 
not all of the possible order 0(e 2 ) terms are retained. Since the HBE model shares the same 
conservation laws and Euler-Poincare structure as EB and PE, and differs from them only 
at order 0(e 2 ), it may be a valid approximation for times longer than the expected order 
0(1 /e) for BE. 



The Kelvin-Noether Theorem. The HBE model ( 4.37 ) possesses the following Kelvin- 
Noether circulation theorem, 

^- I (R + eu R ) ■ dx 3 = - (f bdz, (4.38) 

for any closed curve 74(113) moving with the fluid velocity 113. We compare this result with 
the Kelvin-Noether circulation theorem for PE in equation ( |4.33| ), rewritten as 



— j> (R + eu) -dx 3 = j> (R + eu R + e 2 u D ) ■ dx 3 = - j> bdz. (4.39) 
7t(U3) PE 7f(U3) HBE ZERO 



u D ■ dx 3 = f d\ = 



the e 2 term vanishes, and so the HBE circulation integral differs from that of PE only 
through the differences in buoyancy between the two theories. 

The conservation of potential vorticity on fluid parcels for the HBE model is given by, 



cf. equation ( |4.8| ) 

dl?HBE 

dt 



u 3 • V 3 (?hbe = , where g H BE = V 3 b ■ V 3 x (eu R + R) . (4.40) 
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Combining this with advection of b and the tangential boundary coundition on u 3 yields an 
infinity of conserved quantities, 



for any function These are the Casimir functions for the Lie-Poisson Hamiltonian for- 
mulation of the HBE given in Holm [1996]. 

HBE Discussion. By their construction as Euler-Poincare equations from a Lagrangian 
which possesses the classic fluid symmetries, the HBE conserve integrated energy and con- 
serve potential vorticity on fluid parcels. Their Lie-Poisson Hamiltonian structure endows 
the HBE with the same type of self-consistency that the PE possess (for the same Hamilto- 
nian reason). After all, the conservation laws in both HBE and PE are not accidental. They 
correspond to symmetries of the Hamiltonian or Lagrangian for the fluid motion under con- 
tinuous group transformations, in accordance with Noether's theorem. In particular, energy 
is conserved because the Hamiltonian in both theories does not depend on time explicitly, 
and potential vorticity is conserved on fluid parcels because the corresponding Hamiltonian 
or Lagrangian is right invariant under the infinite set of transformations that relabel the 
fluid parcels without changing the Eulerian velocity and buoyancy. See, e.g., Salmon [1988] 
for a review of these ideas in the GFD context, as well as Holm, Marsden and Ratiu [1998a, b] 
and the earlier sections of the present paper for the general context for such results. 

The vector fields which generate these relabeling transformations turn out to be the 
steady flows of the HBE and PE models. By definition, these steady flows leave invariant 
the Eulerian velocity and buoyancy as they move the Lagrangian fluid parcels along the 
flow lines. Hence, as a direct consequence of their shared Hamiltonian structure, the steady 
flows of both HBE and PE are relative equilibria. That is, steady HBE and PE flows are 
critical points of a sum of conserved quantities, including the (constrained) Hamiltonian. 
This shared critical-point property enables one, for example, to use the Lyapunov method 
to investigate stability of relative equilibrium solutions of HBE and PE. See Holm and 
Long [1989] for an application of the Lyapunov method in the Hamiltonian framework to 
the stability of PE relative equilibria. According to the Lyapunov method, convexity of the 
constrained Hamiltonian at its critical point (the relative equilibrium) is sufficient to provide 
a norm that limits the departure of the solution from equilibrium under perturbations. See, 
e.g., Abarbanel et al. [1986] for applications of this method to the Euler equations for 
incompressible fluid dynamics and Holm et al. [1985] for a range of other applications in 
fluid and plasma theories. 

Thus, the HBE arise as Euler-Poincare equations and possess the same Lie-Poisson 
Hamiltonian structure as EB and PE, and differ in their Hamiltonian and conservation laws 
by small terms of order 0(e 2 ). Moreover, the HBE conservation laws are fundamentally of 
the same nature as those of the EB equations and the PE from which they descend. These 
conserved quantities — particularly the quadratic conserved quantities — may eventually 
be useful measures of the deviations of the HBE solutions from EB and PE solutions under 
time evolution starting from identical initial conditions. 

4.5 Remarks on Two-dimensional Fluid Models in GFD 

The search for simpler dynamics than those of the primitive equations natually leads to 
considerations of two-dimensional fluid models. This certainly holds for applications in GFD, 
because the aspect ratio of the domain (a) and the Rossby number (e) of the rotating flow 
are often small in these applications. Many treatments of two-dimensional GFD models have 
been given using asymptotic expansion methods in Hamilton's principle, see, e.g., Salmon 




(4.41) 
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[1983, 1985, 1988]. These treatments tend to focus especially on the rotating shallow water 
(RSW) equations, their quasigeostrophic (QG) approximation, and certain intermediate 
approximations, such as the semigeostrophic (SG) equations (Eliassen [1949], Hoskins [1975], 
Cullen and Purser [1989], Holm, Lifschitz and Norbury [1998]) and the Salmon [1985] L x 
model. A unified derivation of the RSW, L\, QG and SG equations using Hamilton's 
principle asymptotics is given in Allen and Holm [1996]. This paper also derives as Euler- 
Poincare equations a new class of "extended geostrophic" (EG) models. The EG models 
include nonlocal corrections to the ageostrophic velocity which could produce more accurate 
models than the L\, QG and SG approximations of the RSW equations. 

There are also three dimensional versions of the QG and SG equations, and recently a 
continuously stratified L\ model was derived in Allen, Holm and Newberger [1998] through 
the use of Hamilton's principle asymptotics and the Euler-Poincare theory. For the suite of 
idealized, oceanographic, moderate Rossby number, mesoscale flow test problems in Allen 
and Newberger [1993], this continuously stratified L\ model produces generally accurate ap- 
proximate solutions. These solutions are not quite as accurate as those from the BEM/HBE 
or BE models, but are substantially more accurate than those from three dimensional SG 
or QG. 

Due to their wide applicability in GFD, the properties of the two dimensional QG equa- 
tions have been studied extensively. Wcinstcin [1983] wrote down a Lie-Poisson bracket for 
them in preparation for studying stability of quasigeostrophic equilibria. The Hamiltonian 
structure and nonlinear stability of the equilibrium solutions for the QG system and its 
variants has been thoroughly explored. For references, see Marsden and Weinstein [1982], 
Weinstein [1983] and Holm et al. [1985]. See also the introduction and bibliography of 
Marsden et al. [1983] for a guide to some of the history and literature of this subject. A dis- 
cussion of the geodesic properties of the QG equations in the framework of Eulcr-Poincarc 
theory is given in Holm, Marsden and Ratiu [1998a, b]. A related discussion of QG in both 
two and three dimensions is given from the viewpoint of Hamilton's principle asymptotics 
in Holm and Zcitlin [1997]. 

Formulae showing the asymptotic expansion relationships among the Lagrangians for 
the various GFD models are summarized in Tables (4.1) and (4.2). In the next two sections, 
we turn our attention to dealing with the mean effects of rapid fluctuations in GFD. 
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Table 4.1: Table 1. Successive GFD approximations in Hamilton's principle. 



feuicr = J d 3 x D(l + b)(^B^ji + ^\u\ 2 + ^a 2 w 2 ^j 

Rotation KincticWgy 



v v ' Constraint 



Potential Energy 

^Euior —* Ieb, for small buoyancy, b = 0(e). 
Ieb —* Ipe, for small aspect ratio, a 2 = 0(e). 

Ipe — ► I hbe . for horizontal velocity decomposition, u = z x V?/> + eVx = ur + cud, 
and |u| 2 -> |ur| 2 in ? PE . 

Z hbe — » Zeg, upon rearranging KE in / hbe and decomposing horizontal velocity as 

B E M B E M 

u = Ui + eu 2 , where Ui=zx V<^ with 

/*° 

<f>(x3,t) = <f> s (x,y,t) + / cfe'&, 

•/ z 

i.e., d(j)/dz = — b and where u 2 is the prescribed function, 

u 2 = r(ui • V)z x ui - arV^-V 2 )- 1 ^,^)) -/fr/mi, 

with ^> = /i — &i + V 2 </>. The constants r, a, /?, and 7 are free parameters and the 
functions /1 and b\ denote prescribed order 0(e) Coriolis and topography variations. 

'eg —* h, for horizontal velocity decomposition, u = Ui = z x and dropping terms 
of order 0(e 2 ) in /eg- 

h — * /qg, on dropping terms of order 0(e 2 ) in the Euler-Poincare equations for l\. 
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Table 4.2: Table 2. Nondimensional Euler Poincare Lagrangians at successive 
levels of approximation via asymptotic expansions. 



^Eulor 



D(l + b) R(x) ■ u + - |u| 2 + -a\S --)-p(D-l) 



'EB 



PE = / 



D R • u + -|u| 2 + -a 2 w 2 -bz\ -p(D- 1) 



D R-u+-|u| 2 -6z ] - p(D-l) 



/ hbe = / rf 3 a; 



^eg = / rf 3 a; 



L> R-u+-|u-eu D | 2 -6zj -p(D-l) 



£>( (R + eui +e 2 u 2 ) • u -| Ul + 7 eu 2 | 2 - bz 1 - p(D - 1) 



ii = / d 3 x 



D (R + eu 1 )-u--|u 1 | 2 -6z -p(D-l) 



^qg/aw = / d 2 ^/ rfz 



D I R • u + -u • (1 - C(z) A-^u I - p(D - 1) 



where 



and £? = for QG. 



2G 



5 Generalized Lagrangian Mean (GLM) Equations 



The GLM theory of Andrews and Mclntyre [1978a] is a hybrid Eulerian-Lagrangian de- 
scription in which Langrangian-mean flow quantities satisfy equations expressed in Eulcrian 
form. A related set of equations was introduced by Craik and Leibovich [1976] in their 
study of Langmuir vortices driven by a prescribed surface wave field. The GLM equations 
are extended from prescribed fluctuation properties to a theory of self-consistent Hamilto- 
nian dynamics of wave, mean-flow interaction for a rotating stratified incompressible fluid 
in Gjaja and Holm [1996]. 



GLM Approximations. In GLM theory, one decomposes the fluid trajectory at fixed 
Lagrangian label £ A , A = 1, 2, 3, as follows, 



x^(£ A ,t) = x(e£ A ,et) + a£(x,i), with £ = and x-£ = 0. 



(5.1) 



where scaling with e denotes slow Lagrangian space and time dependence, and an overbar 
denotes an appropriate time average at fixed Eulerian position. For example, overbar may 
denote the average over the rapid oscillation phase of a single-frequency wave displacement 
£ of amplitude a relative to its wavelength. Thus, the displacement £(x, t) associated 
with such waves has zero Eulerian mean, £(x, t) = 0. Superscript £ on a function denotes 
its evaluation at the displaced Eulerian position associated with the rapidly fluctuating 
component of the fluid parcel displacement. Thus, e.g., x^ = x + a£ and u^(x) = u(x^), for 
a function u. 

The GLM operator ( ) averages over parcels at the displaced positions x^ = x + a£ 
and produces slow Eulerian space and time dependence. This defines the Lagrangian mean 
velocity: 



u L (ex, et) = u(x, ff = u? (x, t) = u(x«, t) , 



(5.2) 



where unadorned overbar denotes the Eulerian average and scaling with e denotes slow 
dependence. Thus, the GLM description associates to an instantaneous Eulerian velocity 
field u(x£,t) a unique Lagrangian mean velocity, written (with a slight abuse of notation) 
as u L (ex, et), such that when a fluid parcel at x^ = x + a£ moves at its velocity u(x^, t), a 
fictional parcel at x is moving at velocity u L (ex, et). Hence, 



u(x + a£,t) = (— +u L • V) [x + o£(x,t) 



(5.3) 



where D L = d/dt\ x + u L • V is the material derivative with respect to the slowly varying 
Lagrangian mean velocity, u L (ex, et). 

In GLM theory one finds the following basic Eulerian relations and definitions, 



u ? (x,i) = u(x ? ,t) 
u L (ex,et) 



u L (ex, et) + au'(x,t), 
u(x + a£,t), 



(5.4) 



R«(x) = R(x«) = R L (ex) + aR'(£), 



R L (ex) ee R(x + a^), 

i7= o =rJ, 
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where R denotes the vector potential for the Coriolis parameter, as before. The basic 
identity used in deriving these formulae is 



(I +" f -5f)'<**>-( 



U 



d_ 
dx 



f(x + a£,t), 



(5.5) 



which holds for any differentiable function /. This identity may be shown by taking the 
partial time derivative at constant £ A of the decomposition ( |5.1| ) and using the chain rule, 
cf. Andrews and Mclntyre [1978a]. 



The Lagrangian. We return to Hamilton's principle with Lagrangian ( |1.26| ) for the Euler- 
Boussinesq equations. This is expressed in terms of unsealed (i.e., dimensional) Eulerian 
instantaneous quantities in the GLM notation as 



6 



EB — 



dt 



d 3 3 



|u«(x,i)| 2 -b*gz+ u«(x,<) • R«(x) 



(5.6) 



where g is the constant acceleration of gravity and 

£>* (x, t) = D(x + a£, t) = dot (#j + a^j) , (5.7) 
which is a cubic expression in a. We denote the corresponding pressure decomposition as, 



(x, t) — p(x + a£,t) = p L (ex, et) + op hj (ex, et)pj (x, t) . 



(5.8) 



After expanding D% in powers of a, Gjaja and Holm [1996] average over the rapid space 
and time scales in the action (^^) for the Euler-Boussinesq equations while keeping the 
Lagrangian coordinates l A and et fixed. (This is the Lagrangian mean of the action.) 

Remark. Note that averaging in this setting is a formal operation associated with the 
addition of a new degree of freedom describing the rapid fluctuations. Thus, averaging in 
itself does not entail any approximations. The approximations occur next, in the truncations 
of the expansions of the averaged action in powers of the small parameters e and a. 



The GLM dynamics follows upon making the decompositions in (|5.l|)-(|5T5|) and (5.7) 



(5.8), averaging the action (5.6) and assuming that the rapidly fluctuating displacement £ is 
a prescribed function of x and t, which satisfies the transversality condition given in (5.1C) 
below. 

The averaged, truncated action is (taking R(x) = l~i x x with constant rotation frequency 
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f2 for simplicity) 

&GLM = 



dt / d 6 x( D 




)(^(x,t)) 



gz 



u L ■ (n x x) + « 2 (§ + (n L • v)$) • (n x 



I ' • — ("7^ - 



2 <9x l V da;-? 



dxi ) 



a hApi 



(5.9) 



+0(a 4 ).| 

Note that the buoyancy b is a function of the Lagrangian labels i A which is held hxed during 
the averaging. 



The variation of &glm = J dtLoLU m equation (5.9) with respect to h\ at fixed x and 
t yields the transversality condition, 



(5.10) 



When £ and p\ are single-frequency wave oscillations, this condition implies that the wave 
amplitude is transverse to the wave vector; hence the name "transversality condition." 
Gjaja and Holm [1996] show that this condition is required for the Euler-Poincare equations 
resulting from averaging in Hamilton's principle to be consistent with applying the method 
of averages directly to the Euler-Boussinesq equations. 

Next, the variation of &glm with respect to p L at fixed x and t gives 



„ , a 2 d f..d& , -d^\ 
D = H ( f l — (o —2- 

2 dx l V 8x3 6 dxi J ' 



(5.11) 



where the second term on the right side is order 0(a 2 e), and thus is negligible at order 
0(a 2 ). This follows because mean quantities only depend on slow space and slow time, for 
which d/dx 3 = ed/d(ex J ). 



The Euler Poincare Equations. The Euler-Poincare equation for the action &glm 
results in the GLM motion equation, 

d 



m 



det U X I dex X D 



1 



bgz 



(5.12) 



+ 



d 
9ex 



p L + \u L \ 2 



1- 



-rrL 



- ^u? ■ u£ — u£ • (n X £) — a 2 p • u 



= 0, 



where the Lagrangian mean momentum m is defined as 



m ee 



Su 1 



= D [u L + (O x x)] - a 2 p, 



(5.13) 



and (leaving a 2 explicit) p is the "pseudomomentum density," defined by, cf. Andrews and 
Mclntyre [1978a], 

'dl A " 



p = - D (u) + i? )Vf? and D = det 



dx 3 



(5.14) 



Equations (5.11) and (5.12) (and implicitly ( 5.1 0| ) ) are the equations of the GLM theory for 
incompressible flow discussed in Andrews and Mclntyre [1978a]. 
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The Kelvin— Noether Theorem. The Kelvin-Noether circulation theorem for the GLM 
motion equation (5.12) states that 

— f m-dx^--/ bdz, (5.15) 

where the curve ^f(et) moves with the Lagrangian mean fluid velocity u L . By Stokes' theo- 
rem, the GLM equations generate circulation of m around j(et) whenever the gradient of 
bouyancy is not vertical. The conservation of potential vorticity on fluid parcels for the 
GLM equations is given by 



D qcLM = , where q GL 



1 db / d m 
w = Dd^i ' ( dac X iJ 



(5.16) 



The constraint D = 1 is imposed at order 0(a 2 ) in Hamilton's principle for the action 
(^9) by the slow component of the pressure. This condition holds at order 0(a 2 ) when 
V ex • u L = 0, since D satisfies dD/det + V ex • Du L = 0. 

Alternative Derivation of the Kelvin Noet her Circulation Theorem for GLM. 

The unapproximated Euler-Boussinesq equations ( 4.27 ) for a rotating stratified incompress- 
ible fluid are, in the GLM notation, 



d 



<9x« 



u« u« x 2J7(x«) 

d 



(5.17) 



6« = 0, ^ • u« - 0, 20(x« ) 



<9x« 



R(xf). 



Being Euler-Poincare, these equations admit the following Kelvin-Noether circulation the- 
orem, cf. equation (4.29), 

d 



^ I (u(x*,t)+R(x*)) -d^ = -gl 6(x«,t)dz« 

for a contour j(t) which moves with the fluid. 

The time average of the Kelvin circulation integral is expressed as 



(5.18) 



I(et) 



I (u(x«,i) +R(x«)) • dx« 

J 1 (t) v ' 



7(i) 



(u L + R + au l + aR l ) ■ (dx + a d£) 



7(et) L 



(u + R ) • dx + a 2 (u l + R') ■ d£ 



(5.19) 



where the contour j(et) moves wi th ve locity u £ , since it follows the fluid parcels as the 
average is taken. Using equations ( 5.1 3| ) and ( 5.14| ), we rewrite I(et) as 



I(et)=f (u i + R i -a 2 p/D) • dx = f m • dx 



(5.20) 



which re-introduces Lagrangian mean momentum m(ex, et) and the pseudomomentum den- 
sity of the rapid motion, p(ex, et). In terms of these quantities, we may write the Lagrangian 
mean of the Kelvin-Noether circulation theorem for the Euler-Boussinesq equations as 



dj, d 
dt det 



m ■ dx 



7(et) 



g <2) 6(x£, t) dzt 

hit) 



g (p bdz . 



(5.21) 
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This result recovers the Kelvin-Noether circulation theorem (5.15) for the GLM equations 
(5.12) which were derived above as Euler-Poincare equations for the averaged Lagrangian 
&GLM in equation (5.9). 



Geometry of the Sto kes Mean Drift Velocity and the Pseudomomentum. The 

quantity p in equation ( 5.14 ) is an additional slowly varying fluid degree of freedom which 
emerges in the process of averaging to describe the rectified effects of the rapidly varying 
component of the fluid motion acting on the slowly varying component. The rectified effects 
of wave fluctuations in fluids are traditionally discussed in terms of another quantity, namely, 
the Stokes mean drift velocity, defined as 



u 5 = u L - u = V)u' + 0(a 4 ), 



(5.22) 



We define the analogous Stokes mean drift quantity corresponding to the rotation vector 
potential, R, as 



R S = R L — R = a 2 (£ ■ V)R* + 0(a 4 ) , 
The time-averaged contour integral I(et) in equation (|5.19|) may be rewritten as 



(5.23) 



I(et) 



7(et) 



7(et) 



7(et) 



(u + R) • dx + (u + R ) • dx - p • dx 



(u + R) • dx 
(u + R) • dx 



£■ V(u'+R') -dx + a 2 (u l -d£ 



a 2 [ £^[(u l +n l ) ■ dx 



(5.24) 



Thus, the last two terms combine into the time-averaged Lie derivative with respect to £ of 
the total fluctuation circulation (u ; + R z ) • dx. This is a remarkable formula which shows 
the geometric roles of the Stokes mean drift velocity and the pseudomomentum: The Stokes 
mean drift velocity derives from transport, while the pseudomomentum is caused by line- 
element stretching by the fluctuations. For divergence-free fluctuations, we have V • £ = 0, 
and the Stokes mean drift velocity in the rotating frame is a higher order term. In particular, 



a 2 u s 



dt 







dx 1 



dt 



, d 
dtx 1 



^ dt 



(5.25) 



The last step in this calculation follows because mean quantities only depend on slow space 
and slow time. Thus, the Stokes mean drift velocity (i.e., the difference between the La- 
grangian and Eulerian mean velocities) is order 0(e) smaller than the pseudomomentum, 
for divergence- free fluctuations. 



6 Nonlinear Dispersive Modifications of the 
EB Equations and PE 

In generalizing earlier work by Camassa and Holm [1993] from one dimension to n dimen- 
sions, Holm, Marsden and Ratiu [1998a, b] used the Euler-Poincare framework to formulate 
a modified set of Euler equations for ideal homogeneous incompressible fluids. This modi- 
fication introduces nonlinear dispersion into Euler's equations, which is designed physically 
to model nondissipative unresolved rapid fluctuations. This nonlinear dispersion is founded 
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mathematically on the geometrical property that solutions of the basic Euler-Poincare equa- 
tions are geodesies on an underlying group when their Lagrangian is a metric. Here we use 
the Euler-Poincare reduction theorems ( [2.l| ) and (3.1) including advected parameters and 



Hamilton's principle asymptotics in the EB Lagrangian in GLM notation (5.6) to formu- 
late a modified set of Euler-Boussinesq equations that includes nonlinear dispersion along 
with stratification and rotation. In this new modification of the Euler-Boussinesq equations, 
nonlinear dispersion adaptively filters high wavenumbers and thereby enhances stability and 
regularity without compromising either low wavenumber behavior, or geophysical balances. 
Here, 'high' and 'low' refer repectively to wavenumbers greater, or less than the inverse of 
a fundamental length scale a, which parameterizes the nonlinearly dispersion in the model. 
We also present the corresponding nonlinear dispersive modification of the primitive equa- 
tions. We leave it as an open question, whether our nonlinearly dispersive primitive equation 
model will have a slow manifold when dissipation and forcing are included. 

6.1 Higher Dimensional Camassa— Holm Equation. 

The Lagrangian and Action Functionals. As shown in Holm, Marsden and Ratiu 
[1998a, b], the Camassa-Holm (CH) equation (Camassa and Holm [1993]) in n dimensions 
describes geodesic motion on the diffcomorphism group of R" with respect to the metric 
given by the H 1 norm of the Eulerian fluid velocity. This Euler-Poincare equation follows 
from the Lagrangian leu given by the H 1 norm of the fluid velocity u in n dimensions, 
subject to volume preservation (for n 1), namely, 



©ch = j dtl C H = J dtjd n x j|(|u| 2 + a 2 |Vu| 2 ) -p(£>-l)J, 



(6.1) 



in which the parameter a has dimensions of length. We denote (Vu)j = ttj = du % /dx 3 , 
Vu| 2 = u' jV,f = tr(Vu • Vu T ), tr is the trace operation for matrices and the superscript 
(•) T denotes transpose. 

The action ©ch in ( |6 . 1[ ) is the order 0(a 2 ) expression for the mean of the EB action in 
equation ( |5.6| ) in the absence of stratification and rotation, namely 

©ch « J dt J d n x j#i|u«(x,i)| 2 -p(D-l) J, (6.2) 

in which we neglect the order 0(a 2 e) pressure and density corrections discussed in Gjaja 
and Holm [1996] and approximate the mean kinetic energy in a Taylor expansion as follows, 

\WM? = i|u(x + ^,t)| 2 = i|u(x,t) + a£.Vu| 2 + 0(a 4 ) 

= i|u| 2 + ^[TVuP + 0(« 4 ) a I| u | 2 + ^|Vu| 2 . (6.3) 

In the last step, we have also dropped terms of order 0(a 4 ) and assumed isotropy of the 
rapid fluctuations, so that » 8 tJ . As in GLM theory, the length-scale parameter a 
represents the amplitude of the rapidly fluctuating component of the fluid parcel trajectory 
(or its amplitude to length-scale ratio in a nondimcnsional formulation). 
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The Euler Poincare Equations. Varying the action ©ch in (6.1) at fixed x and t gives 



se 



CH 



dt I d n x 

M 



|Vu| : 



p) 6D 



+ (Du - a 2 (divDgrad)u) • Su - (D - l)8p 
+ a 2 [dt(f d n ~ l x (Dn- Vu-Su), 

J JdM 



whose natural boundary conditions on dM are 

u • h = and (h • V)u || h, 



(6.4) 



(6.5) 



where || denotes "parallel to" in the second boundary condition, which of course is not 
imposed when a 2 is absent. (Recall that Su in equation (|6.4[) is arbitrary except for being 



tangent on the boundary. This tangency, along with the second condition in equation (6.5) 
is sufficient for the boundary integral in equation ([O]) to vanish.) 

Unfortunately, the space of divergence free vector fields with the boundary conditions 



(S.5) do not form a Lie algebra unless the boundaries are flat. There are two ways to deal 
with this problem. One is to use the no-slip boundary conditions u = (and hence Su — on 
dM.). This is especially interesting in conjunction with the addition of viscosity. Another is 
to modify these boundary conditions with the addition of a second fundamental form term 
to the right hand side. The latter idea is explored in detail in Holm, Kouranbaeva, Marsden, 
Ratiu, and Shkol lcr [1 998] . In what follows, we imagine that one of these options is chosen. 
By equation ( [f.ll|) , the Euler-Poincare equation for the action ©ch in equation (5.1) is 



d 

— + u • V I v + VjVu J 
dt 



V(y-i|u| a -y|Vuf 



= 0, 



(6.6) 



where 



u - a Au 



1 61 



CH 



D Su 



D=l 



In writing this equation, we have substituted the constraint D = 1, which as before implies 
incompressibility via the continuity equation for D. Requiring the motion equation (6.£) to 
preserve div u = = div v implies a Poisson equation for the pressure p with a Neumann 
boundary condition, which is obtained just as usual in the case of incompressible ideal 
fluid dynamics, by taking the normal component of the motion equation evaluated at the 
boundary. Of course, the n-dimensional extension of the CH equation (6J3) reduces to 
Euler's equation when we set a = 0. The properties of the n-dimensional CH equation (6.6) 
are summarized in Holm, Marsden and Ratiu [1998a, b] and Holm, Kouranbaeva, Marsden, 
Ratiu, and Shkoller [1998] for the ideal case. In particular, equation ( |6.6| ) is shown to be 
the geodesic spray equation for geodesic motion on the group Diff(X>) with respect to the 
metric given by the H 1 norm of the fluid velocity u in n dimensions. See also Chen, Foias, 
Holm, Olson, Titi and Wynne [1998] and Foias, Holm and Titi [1998] for discussions of the 
corresponding viscous, forced case and its connection to mean turbulence closure models. 



Discussion of the CH Equation. The essential idea of the CH equation is that its 
specific momentum (i.e., its momentum per unit mass) is transported by a velocity which 
is smoothed by inverting the elliptic Hclmholtz operator (1 — a 2 A), where a corresponds 
to the length scale at which this smoothing becomes important, i.e., when it becomes of 
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order 0(1). When the smoothing operator (1 — a 2 A) _1 is applied to the transport velocity 
in Euler's equation to produce the CH equation, its effect on length scales smaller than a 
is that steep gradients of the specific momentum v tend not to steepen much further, and 
thin vortex tubes tend not to get much thinner as they are transported. And, its effect 
on length scales that are considerably larger than a is negligible. Hence, the transport of 
vorticity in the CH equation is intermediate between that for the Euler equations in two and 
three dimensions. As for Euler vorticity, the curl of the CH specific momentum is convected 
as an active two form, but its transport velocity is the smoothed, or filtered CH specific 
momentum. 

The effects of this smoothing or filtering of the transport velocity in the CH equation 
can be seen quite clearly from its Fourier spectral representation in the periodic case. In 
this case, we define as the k-th Fourier mode of the specific momentum v = (1 — a 2 A)u 
for the CH equation; so that vj. = (1 + a 2 |k| 2 )uj <; . Then the Fourier spectral representation 
of the CH equation for a periodic three-dimensional domain is expressed as 



x (n x v n ) = 0, (6.7) 



where n is the Leray projection onto Fourier modes transverse to k. (As usual, the Leray 
projection ensures incompressibility. ) In this Fourier spectral representation of the CH 
equation, one sees that the coupling to high modes is suppressed by the denominator when 
1 + a 2 |p| 2 3> 1. Consequently, when |p| > 0(1 /a), the smoothing of the transport veloc- 
ity suppresses the Fourier-mode interaction coefficients. In fact, the CH smoothing of the 
transport velocity suppresses both the forward and backward cascades for wave numbers 
|p| > 0(1/ a), but leaves the Euler dynamics essentially unchanged for smaller wave num- 
bers. The result is that the vortex stretching term in the dynamics of q = curlv is mollified 
in the CH model and so the vortices at high wave numbers will tend to be "shorter and 
fatter" than in the corresponding Euler case for the same initial conditions. 




The Kelvin— Noether Theorem. Since the ri-dimcnsional CH equation (6.6) is Euler 



Poincare, it also has a corresponding Kelvin-Noether circulation theorem. Namely, 

^ = where I(t) = I (u - a 2 Au) • dx = I v-dx, (6.8) 



for any closed curve jt that moves with the fluid velocity u. This expression for the Kelvin- 
Noether circulation of the CH equation in three dimensions is reminiscent of the correspond- 
ing expression ( 5.20| ) in GLM theory involving the pseudomomentum for wave, mean- flow 



interaction theory. This correspondence confirms the physical interpretation of the a 2 term 
in the Kelvin-Noether circulation integral as a Lagrangian mean closure relation for the 
pseudomomentum of the high frequency (i.e., rapidly fluctuating, turbulent) components of 
the flow. In this interpretation, a for the CH equation corresponds to both the amplitude 
of these high frequency components and the length scale at which they become important. 



Energy Conservation for the CH Equation. Legendre transforming the action (6.1) 
gives the following conserved Hamiltonian (still expressed in terms of the velocity, instead 
of the momentum density m = Sl/Su), 



H= d n x 
Jm 



|(|u| 2 + a 2 |Vu| 2 )+K£>-l)]. (6.9) 
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Thus, when evaluated on the constraint manifold D = 1, the Lagrangian and the Hamilto- 
nian for the CH equation coincide in n dimensions. (This, of course, is not unexpected for 
a stationary principle giving rise to geodesic motion.) 

The curl of the 3D Camassa-Holm motion equation ([T(]) yields 

d 

— q = q • Vu — u ■ Vq = [u, q], where q = curl(u — a 2 Au) , (6.10) 

and we have used incompressibility and commutativity of the divergence and Laplacian 
operators. Thus, u is the transport velocity for the generalized vorticity q and the "vortex 
stretching" term q • Vu involves Vu, whose L 2 norm is controlled by the conservation of 
energy in equation ( |6.9[ ). Boundedness of this norm will be useful in future analytical studies 
of the 3D Camassa-Holm equation; for example, in the investigation of the Liapunov stability 
properties of its equilibrium solutions. 



The Riemannian CH Equation. One can formulate the CH equation on a general 
Riemannian manifold, possibly with boundary. Although this formulation will be the subject 
of future papers, we comment on it here because of the importance of spherical geometry, 
in particular, for GFD models. 

6.2 The Euler-Boussinesq a Model (EBa) 

We introduce nonlinear dispersion into the Euler-Boussinesq equations using the Euler- 
Poincare framework and following the example of the CH equations. 



The Lagrangian. To carry this out, we modify the EB Lagrangian ( 4.26 ) by simply 
adding the a 2 term (while dropping other scale factors e, <r, as well as primes and subscripts) 



Seb« = dt Zebq 



dtl d n x 

M 



(6.11) 



/)( i|u| 2 + ^|Vu| 2 + u • R(x) - gbz ) i>i D I i 



This action is the order 0(a 2 ) ap prox imation of the mean of the EB action in equation 
(5.6), using the Taylor expansion (6^3) and neglecting order 0(a 2 e) pressure and density 
corrections. (An order 0(a 2 ) term involving (£ • Vu)(£ • VR) is also dropped, since it makes 
only a negligible contribution in the resulting motion equation.) 



The Euler Poincare Equations. Varying this action at fixed x and t gives 
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EBa 



dt I d n x 

M 



1 



| Vu| 2 + u • R(x) - gbz - p) 8D 



2' 1 2 

DgzSb -(D- l)5p + (Du - a 2 (divDgrad)u) • Su 
+ a 2 / dt I d n ^x (Dn ■ Vu • Su) , 

J JdM 



(6.12) 



where the natural boundary conditions are again given in (6.5). The corresponding Eulcr- 
Poincare equation for the action Sebq in equation ( 6.11 ) is, 



0_ 

di 



u • V ) v — u x curlR + vjVu^ + gbz 

+^ ! p-^|u| 2 -^|Vu| 2 ) =0, 



(6.13) 
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where 



u - a Au, 



db _ Q ±_ fd_ 

dt ~ ' di ~~ \dt 



V • u = 0, 
u- V 



(6.14) 



Again, we have substituted the constraint _D = 1, which implies incompressibility via the 
continuity equation for D. Relative to the usual EB equations (4.27), the EBa equation 



(5.13) has a smoothed, or filtered transport velocity, since u = (1 — a 2 A) 1 v 



The Kelvin— Noether Theorem. The Kelvin-Noether circulation theorem for the EBa 
equation ( |3.13| ) is, 

d_ 

dt 



(p (v + R) • dx = — &> bdz 

•'7t( u ) "'7t( u ) 



(6.15) 



where the curve j t (u) moves with the fluid velocity u. (The two Kelvin theorems in equations 
( |6.15 ) and (4.29) differ in their definitions of v.) By Stokes' theorem, the EBa equations 
generate circulation of v + R around 7t(u) whenever the gradient of bouyancy is not vertical. 
The conservation of potential vorticity on fluid parcels for the EBa equations is given by 



/EBa 



at 



u • V qEBa = , where <?EBa = Vfe • V x (v + R) . 



The curl of the EBa motion equation gives, 
d 



dt 



u • Vq + q • Vu + gVb x z, where q = curl(u — a 2 Au + R) 



(6.16) 



(6.17) 



This the usual expression for transport, stretching and creation of vorticity in a buoyant 
flow, except that here the vortex stretching coefficient is Vu, which is moderated relative to 
the usual EB equations. 



Energy Conservation. The EBa equations (6.13) conserve the following Hamiltonian 
(found, e.g., by Legendre transforming the Lagrangian Iebo in equation (5.11)). Namely, 



^EBa = 



M 



y|Vu| 



Dgbz + p(D - 1) 



The corresponding conserved energy is 



E 



EBc 



Hi 



EBc 



15=1 



d n x 



ri 



M 



a. 

Y 



Vu| : 



gbz 



(6.18) 



(6.19) 



Since the (finite) value of this conserved energy for the EBa model is determined by its initial 
conditions and b is advected (so it has a maximum value in L°°) there is H 1 control of the 
velocity u. This is the effect of the "filtering" of the solution produced by the nonlinear 
dispersion for q / 0. This filtering moderates the growth of instabilities at wavenumbers 
k| > 1/a. So if the EBa model is used as a large eddy simulation (LES) model, one would 
choose the value of a to determine the size of the minimum resolved length scale. The 
filtering by the a term also allows nonlinear Liapunov stability conditions to be formulated 
for equilibrium solutions of the EBa model. This stability result is clear from the work 
of Abarbanel et al. [1986], who introduced the notion of "conditional" Liapunov stability 
for the EB model, using wavenumber conditions that now turn out to be satisfied for the 



EBa model. The Euler fluid equations (4.23) may also be modified analogously to include 
nonlinear dispersion. However, this case is ignored for now, as we proceed to discuss the 
modified primitive equations. 
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6.3 Primitive Equation a Model (PEa) 

In horizontal and vertical components, with v 1 - = v — v"z, V" = (v-z) and curlR = f(x)z, 



the the EBa equations are expressed in nondimcnsional form as, cf. equations ( [4.2q ) and 



( P-13D , 



e^+fu,VV + /2Xu + V 1 ^ = 0, £( r 2 ^+6+^ = 0, (6.20) 
or at az 

d f 9 i _ i ii 9 \ d& _ _ | , chJI 



where 7r=lp--|u| — |Vu| I , v = u a Au . 

Here, e and a are the Rossby number and aspect ratio, respectively, and a has been scaled 
in units of horizontal length scale, L. The leading order balances are still hydrostatic in 



the vertical, and geostrophic in the horizontal. Setting a = in equation (5.20) removes 
the vertical acceleration and thereby produces the primitive equation a model (PEa). We 
expect that the filtering property for a ^ discussed above for the CH and EBa equations 
should make the PEa model much more regular than the ordinary PE, from the viewpoint 
of gravity wave oscillations. The problematic asspects of gravity waves in atmospheric and 
oceanic numerical simulations and data assimilation using the PE model have often been 
addressed by invoking the concept of an idealized "slow manifold" on which gravity waves 
are absent, see, e.g., Lorenz [1992]. However, the existence of a slow manifold has never 
been proven for the PE model. It is an open question, whether the new PEa model will 
have a slow manifold when dissipation and forcing are included. We shall report on this 
matter elsewhere. 



Final Remarks 

In this paper we have shown how asymptotic expansions in Hamilton's principle for the 
Euler-Poincare equations of geophysical fluid dynamics provide an organizing principle for 
many GFD systems and produce a clear, unified understanding of their Kelvin theorems. 
Hamilton's principle asymptotics in the Euler-Poincare setting thus explains the shared 
properties of these GFD models and provides a unified approach to making additional 
approximations and creating new models. In this setting, we have introduced a new class of 
fluid models, called a models, that possess nonlinear dispersion which smooths the transport 
velocity relative to the circulation velocity in the Euler-Poincare equations. The effect of 
this smoothing is to moderate the vortex stretching process while preserving the Kelvin 
circulation theorem for these equations. The efficacy and utility of the a models are yet to be 
determined, but initial studies of them are promising. We expect that one can also perform 
the sort of asymptotics done here on the group level first (to get an approximating group) and 
then perform Euler-Poincare reduction and arrive at the same conclusions. In other words, 
there should be a general principle for asymptotics in Hamilton's principle, which allows 
passage from one group and its corresponding advected quantities to an approximating one. 
Moreover, this process should commute with Lagrangian and Hamiltonian reduction. 
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